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ABSTRACT 

We have modeled the unusual orbital light curve of V1494 Aquilae (Nova Aquilae 1999 No. 2) and 
found that such an unusual orbital light curve can be reproduced when there exist two-armed, spiral 
shocks on the accretion disk. V1494 Aql is a fast classical nova and found to be an eclipsing system 
with the orbital period of 0.1346138 days in the late phase of the nova outburst. Its orbital light curve 
shows a small bump at orbital phase 0.2, a small dip at 0.3, sometimes a small bump at 0.4, and a 
large bump at 0.6 — 0.7 outside eclipse. Such a double- or triple-wave pattern outside eclipse has never 
been observed even though overall patterns look like some supersoft X-ray sources or eclipsing polars. 
We have calculated orbital light curves including the irradiation effects of the accretion disk and the 
companion by the hot white dwarf. These unusual patterns can be reproduced when we assume two- 
armed spiral shocks on the accretion disk. Especially, triple-wave patterns are naturally obtained. This 
result strongly suggests the existence of two-armed spiral shocks on the accretion disk in the late phase 
of the nova outburst. 

Subject headings: accretion: accretion disks — binaries: close — binaries: eclipsing — novae, 
cataclysmic variables — stars: individual (V1494 Aquilae) 



1. INTRODUCTION 

Angular momentum transport plays an essential role in 
accretion disks of cataclysmic variables. Two mechanisms 
have been proposed so far: one is the turbulent viscosity as 
adopted in the standard accretion disk model of Shakura 
& Sunyaev (1973) and the other is the direct dissipa- 
tion by tidal spiral shocks on the accretion disk as first 
demonstrated by Sawada, Matsuda, & Hachisu (1986). 
The turbulent viscosity is a local physical process while 
the tidal spiral shocks have a global structure on the ac- 
cretion disk. Therefore, we have a chance to observe a 
global spiral shock structure when they play an essential 
role in the angular momentum transport of the accretion 
disk. Such an observational evidence first came from the 
Doppler maps of the dwarf nova IP Peg outburst (e.g., 
Steeghs, Harlaftis, & Horne 1997). We have long believed 
that tidal spiral structures can be detected even in orbital 
light curves of cataclysmic variables if the spiral patterns 
are prominent. At last we find such an evidence of spiral 
patterns on the accretion disk from the orbital light curves 
of V1494 Aquilae. 

V1494 Aquilae (Nova Aquilae 1999 No. 2) was discovered 
by A. Pereira on 1999 December 1.785 UT at my ^ 6.0 
(Pereira et al. 1999). It reached the visual maxi- 
mum brightness of my ^ 4.0 on Dec 3.4 UT (to =JD 
2,451, 515. 9±0.1) and decayed at the rate of = 6.6 ± 0.5 



days and is = 16 ± 0.5 days (Kiss & Thomson 2000). 
Early phase spectra were taken by Fujii (1999) and Ayani 
(1999), which show P-Cyg profiles of hydrogen Balmer 
lines with a blueshifted component of 1020 km s"^ and 
1200 km s"\ respectively. Thus, V1494 Aql has been 
classified as a fast nova. Kato et al. (2004) have sum- 
marized overall development of the V1494 Aql light curve 
until autumn of 2003 (for about 3 years). Detailed spectral 
features have been reported by lijima & Esenoglu (2003). 

A short-term periodic modulation was first reported by 
Novak et al. (2000), who found 0.03 mag variations with 
a period of 0.0627 days from their 2000 June observation. 
Retter et al. (2000) analyzed 31 nights of CCD photome- 
try during 2000 June- August and suggested a periodicity 
of 0.13467(2) days. The full amplitude of the variation in 
the i?-band increased from 0.03 mag in June to 0.07 mag 
in August. Their folded light curve shows a double- wave 
structure with a shallower dip at phase 0.5, with about half 
the amphtude of the main periodicity. Bos et al. (2001) 
reported, based on unfiltered and i?-band CCD photome- 
try obtained on 12 nights during 2001 June- July, a robust 
change in the shape and amplitude of the 0.13467 days 
period. It had an eclipse shape with depth about 0.5 mag. 
A second shallow eclipse (about 0.1 mag deep) at phase 
0.5 can be seen. Barsukova & Goranskii (2003) refined, 
based on V-band photometry during June and September 
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of 2002, the orbital period of 0.1346141(5) days. They 
pointed out that the orbital light curve is quite unusual 
in the sense that it is not similar to ordinary cataclysmic 
variables with a hot spot on the accretion disk. Pavlenko, 
Dudka, & Baklanov (1999) made a multicolor photometry 
and concluded that the eclipse depth is deeper in longer 
wave lengths (deeper in /-band than in F-band). Pavlenko 
et al. (1999) further discussed that such an eclipsing 
characteristic (deeper in longer wave lengths) cannot be 
explained by the accretion disk or the irradiation effect 
of the companion because these two light sources have the 
opposite nature (deeper in short wave lengths). Thus, they 
suggested, as a model of V1494 Aql orbital light curve, a 
self eclipsing accretion column in a magnetic polar. 

Kato et al. (2004) provided a time-development of the 
orbital light curves among 2001 November-December, 2002 
August, and 2003 June. It sometimes shows a triple-wave 
structure as well as a double-wave structure. Kato et al. 
suggested that some structures (probably on the accretion 
disk) fixed in the binary rotational frame are responsi- 
ble for the orbital light curve because of the stability of 
out-of-eclipsc light curve patterns. Very recently. Kiss, 
Csak, & Derekas (2004) pointed out that a close compan- 
ion to V1494 Aql (located 1'.'4 southwest) is brighter than 
V1494 Aql itself in the very late phase of the outburst. 
They corrected the brightness of the orbital light curve 
and found that the depth of the eclipse is about twice 
deeper in magnitudes than before correction. 

The depth of the primary eclipse of V1494 Aql has be- 
come deeper as the nova decayed. Such a feature was 
also observed in the recurrent nova CI Aql 2000 outburst 
(e.g., Matsumoto et al. 2001), in which the irradiation 
of the accretion disk plays an essential role (Hachisu & 
Kato 2003a). Therefore, we expect the accretion disk in 
V1494 Aql is also responsible for the unusual wavy struc- 
ture of the orbital light curve. In this Letter, we model 
the orbital light curve. Almost the same orbital light 
curve model as in the supersoft X-ray sources is adopted 
(e.g., Schandl, Meyer-Hofmeister, & Meyer 1997; Hachisu 
& Kato 2003b, c), which includes the irradiation effects of 
the accretion disk and of the companion. In addition, we 
assume two-armed spiral structures on the accretion disk 
to reproduce a triple-wave structure outside eclipse. In §2, 
our numerical model for V1494 Aql is briefly introduced 
and summarized. The numerical results are given in §3. 
Discussion follows in §4. 

2. THE LIGHT CURVE MODEL 

Our binary model is illustrated in Figure 1, which con- 
sists of a 0.3 Mq main-sequence star (MS) filling its Roche 
lobe, a 1.0 Mq white dwarf (WD), and a disk around the 
WD. The mass of the WD is roughly estimated from the 
decline rate of the decay phase (e.g., Hachisu et al. 2000; 
Hachisu & Kato 2001a, b), details of which will be pub- 
lished elsewhere. The mass of the companion is estimated 
from the mass-period relation for cataclysmic variables 
(e.g., Patterson 1984; Warner 1995). A circular orbit 
is assumed. Its ephemeris has recently been refined by 
Kato et al. (2004). We use this ephemeris, i.e., 

t(BJD) = 2, 452, 458.3230 -|- 0.1346138 x E, (1) 

at eclipse minima. 



We also assume that the surfaces of the WD, the MS 
companion, and the accretion disk emit photons as a black- 
body at a local temperature which varies with position. 
For the basic structure of the accretion disk, we assume 
an axi-symmetric structure with the size and thickness of 
-Rdisk = OiRl , (2) 
and ^ 

/i = /3i?disk (^;^) , (3) 

where -Rdisk is the outer edge of the accretion disk, 
is the effective radius of the inner critical Roche lobe for 
the WD component, h is the height of the surface from 
the equatorial plane, and w is the distance on the equa- 
torial plane from the center of the WD. Here, we adopt 
u = 2. We obtain our modeled two-armed spiral struc- 
tures by multiplying h with -Zheight) as defined below: 



zi = max 1, 



Z2 = max I 1, 



^ (oj/iidisk - exp{-ri{(f) - J)))2 + £2 
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(4) 



(5) 

-Zheight = max(zi, 22). (6) 

We further have a slope at the disk edge defined by 
z - 10{z - 0.25) (cc/ i?disk - 0.9) for 0.9 < tn/i?disk < 1-0. 
The above various disk parameters are assumed to be e = 

0. 1, = ^2 = 0.40, r] = 0.41, S = 65°, o = 0.8, fi = 0.095 
for the disk shape of Figure 1. Here, ^1 and ^2 specify 
the amplitudes of two spirals so that ^1 = ^2 (= 0.40) 
means just an anti-symmetric structure of spirals with the 
same height, rj determines the inverse of the pitch angle 
of logarithmic spirals, 6 is the position angle of the spi- 
rals against the binary components, e denotes the width 
of the spiral pattern and represents the height of the spiral 
against the thickness of the disk together with ^1 and ^2, 

1. e., -Zheight = max(^i, ^2)/e ~ 4 at the edge of the disk for 
the disk shape of Figure 1. In our light curve model, we 
mainly change four parameters, i.e., the disk thickness (/3), 
the inverse of the pitch angle of spirals {rj), the position 
angle (S), and the heights of two spirals (^1 = ^2)- The 
other parameters that specify the disk shape are all fixed 
to be the above values. These parameters are roughly de- 
termined to mimic the results of 3D simulation model of 
accretion disks (e.g. Makita, Miyawaki, & Matsuda 2000). 

The luminosity of the WD is assumed to be 750 — 
3,000 Lq because we do not know the exact WD lumi- 
nosity at the time of observation. The disk and the com- 
panion star are strongly irradiated by the hot WD. The 
surfaces of the WD, the disk, and the companion star are 
divided into many patches as shown in Figure 1. Here 
we assume that each patch emits photons as a blackbody 
(with a single temperature). Each patch of the disk or of 
the companion is irradiated by visible (front side) patches 
of the hot WD. The total luminosity of the irradiated disk 
and companion is calculated by summing up the contri- 
butions from all patches. The irradiation efficiency is the 
same as that of Schandl et al. (1997), i.e., 50% of the irra- 
diated energy is emitted by photon but the residual 50% is 
converted into thermal energy of gas. The accretion lumi- 
nosity of the disk is also numerically included, although its 
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contribution to the optical light is much smaller than that 
of the irradiation effect (see discussion of Hachisu & Kato 
2001a, b; Hachisu, Kato, & Schaefer 2003). The numeri- 
cal method adopted here was fully described in Hachisu 
& Kato (2GGlb, 2GG3a,b,c). The inclination angle (i) of 
the binary is a free parameter that is determined from our 
light curve fitting. The original temperature of the com- 
panion star is assumed to be 3,000 K, but the shape of the 
orbital light curve is hardly changed even if we take 2,000 
K or 4, GOG K for the original temperature of the MS com- 
panion. The adopted system parameters are summarized 
in Table 1. 

3. NUMERICAL RESULTS 

The best fit model is plotted in Figure 2 by a thick solid 
line. The orbital modulations of V1494 Aql have been 
reported by several groups. Here, we adopt data of two 
groups (Kato et al. 2004; Kiss et al. 2004) to fit with 
our modeled light curves. Kato et al.'s (square) data have 
been corrected by using Kiss et al.'s (circle) data, because 
Kato et al.'s original data include the light from the close 
companion. 

We have changed five parameters (/3, rj, 6, = ^2, i) 
independently and calculated the total of ~ 1300 orbital 
light curves. It takes about 6 hours to calculate one orbital 
light curve on a 2.4 GHz Pentium 4 processor, because the 
surface patch elements arc 64 x 128 (61 x (/)) for the com- 
panion surface, 64 x 128 x 2 (tu x <j>x (up and down sides)) 
for the disk, and 16 x 32 (6* x 0) for the WD, and the 
total of 129 steps for one orbital period. We have used 5 
CPUs, so that it took the total computation time of about 
1600 hours, i.e., (6 x 1300)/5 CPUs. We have found, by 
the least square method, the best fit model for out-of- 
eclipse (orbital phase 0.1 - 0.9), i.e., /3 = 0.095, rj = 0.41, 
S = 65°, ^1 = ^2 = 0.40, i = 78.5°. For the best fit 
model, we increase the number of patches for the disk up 
to 128 x 256 X 2 (tu x ^x (up and down sides)) in order to 
calculate a rather smooth light curve. Such a high quality 
light curve is shown in Figure 2. 

The triple-wave structure of out-of-eclipse orbital light 
curve is reasonably and naturally reproduced with our spi- 
ral shock pattern model. Both the thickness of the accre- 
tion disk and the pitch angle of spirals adopted here are 
roughly consistent with the 3D simulations (e.g., Makita 
et al. 2000). The present results strongly support 
the two-armed spiral pattern on the accretion disk. We 
hope a spectroscopic confirmation of this spiral pattern by 
Doppler maps for V1494 Aql. 

4. DISCUSSION 

The orbital light curve varies from night to night as 

shown in Figure 3 of Kato et al. (2004). Such a scatter 
is also clearly shown in the present Figures 2 and 3 for 
Kiss et al.'s (2004) data. To reproduce this variation, we 
have changed parameter ^1 and ^2 independently (see Fig. 
3). The variation appears in orbital phase of 0.1 — 0.3 and 
the upper and lower bounds for these variations can be 
well reproduced by the models of ^1 = 0.26, ^2 = 0.44 for 
the upper bound, and ^1 = 0.44, ^2 = 0.24 for the lower 
bound. Here, ^1 and ^2 represent the height of the spiral 



arm in the rear and front side on the disk of Figure 1, re- 
spectively. The other parameters are the same as those in 
Figures 1 and 2. 

It is remarkable that the peak at orbital phase ~ 0.65 is 
rather stable and never disappears. This is consistent with 
the robustness of the period for the out-of-eclipse shape as 
emphasized by Kato et al. (2004). 

Pavlenko et al. (1999) discussed the main light source 
for the wavy structures of the out-of-eclipse orbital light 
curve. They argued, as a light source, (1) accretion disk, 
(2) ellipsoidal shape of the companion, (3) reflection (irra- 
diation) effect of the companion by the hot white dwarf, 
(4) intermediate polar activity, and (5) polar activity self- 
eclipsed by the accretion column. They rejected possibili- 
ties (l)-(4) and suggested a viable model of polar activity. 
Their reasons rejecting possibilities (2) and (4) are rea- 
sonable but not for possibilities (1) and (3) as discussed 
below. 

Pavlenko et al. (1999) rejected the possibility of accre- 
tion disk by pointing out two reasons: (i) the maximum 
duration of the eclipse of the disk limited by the Roche 
lobe size and it cannot exceed a quarter of the orbital 
phase while the overall duration of the eclipse is as large as 
0.45 of the orbital period, (ii) The peak of light emission 
of the accretion disk lies in the blue region of the spec- 
trum, so that the amplitude of the eclipse should increase 
as the wavelength becomes shorter ( "as the wave length is 
increased" in their text is mistaken). However, they ob- 
served the opposite trend: the amplitude in the red region 
is much more deeper than in the visual region. As for 
the first conjecture (i), based on the accretion disk model 
together with the irradiation effects of the disk and the 
companion, we have already constructed the orbital light 
curves that match well the observational data. Second, we 
should be careful with the above statement (ii) because 
the depth of the primary eclipse varies from night to night 
as seen in Figure 3 of Kato et al. (2004) and its depth 
changes as large as ARc ~ 0.4 between different periods of 
the observations as shown in Figures 2 and 3. Pavlenko et 
al.'s (2003) R and / light curves are not simultaneous ones 
but taken in different periods {R is earlier than /) . The V 
light curve is reconstructed from Barsukova & Goranskii's 
(2003) data. 

They also rejected possibility (3) based on the same 
statement as (ii) discussed above. However, their multi- 
color results may simply suggest the fact that the eclipse 
depth varies from period (night) to period (night) . There- 
fore, we cannot conclude statement (ii) only from the dif- 
ferent depths at the different periods. 

For the assumed WD luminosity of iwD — 3,000 Lq, 
the distance modulus is obtained to be (to — M)ii = 12.05 
and the corresponding distance is estimated to be d 1.4 
kpc (see Table 1). However, this does not mean the real 
distance to V1494 Aql because we do not know the real lu- 
minosity of LwD at the time of the two observations. For 
instance, if we adopt Lwd = 750 Lq, the shape of the or- 
bital light curve hardly changes but the distance modulus 
becomes (to — M)fi — 11.45 and ~ 1 kpc. 

We thank L. Kiss for providing us their V1494 Aql data. 
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Table 1 



Adopted system parameters of V1494 Aquilae 



parameter 


symbol 


present work 


inclination angle 


i 


78.5° 


WD mass 


MwD 


1.0 Mq 


WD luminosity 


^WD 


750 - 3,000 Lq 


MS mass 


Mms 


0.3 Mq 


mass accretion rate 


M 


0.3 X 10"^ Mq yr-i 


MS temperature 


'2\lS.Qrg 


3,000 K 


distance modulus 


(m - M)r 


11.45 - 12.05 


color excess 


E(B - V) 


0.60' 


7?- band absorption 


Ar 


1.39'= 


distance 


d 


1.0 - 1.4 kpc 



^taken from lijima Sz Esenoglu (2003) 

^from extinction law of Rieke Lebofsky (1985) 
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Fig. 1. — Our V1494 Aql model is illustrated. The cool component {far and left side) is an MS companion (0.3 Mq) filling up its inner 
critical Roche lobe. The north and south polar areas of the cool component are irradiated by the hot component (1.0 Mq WD, near and 
right side). The separation is a = 1.21 Rq; the effective radii of the inner critical Roche lobes are R* = 0.59 Rq, and R2 = R2 = 0.34 Rq, 
for the primary WD and the secondary MS companion, respectively. A two-armed spiral pattern is shown on the a<;cretion disk. The WD 
surface is artificially enlarged to easily see it. 



6 



Hachisu et al. 



15 - 



T 1 1 1- 



T 1 1 r 



T 1 1 1- 



VI 494 Aql 1.0 Mq WD + 0.3 Mq MS (m-M)R=12.05 i=78.5" 

Kiss et al.'s data (2003 August) 
□ Kato et al.'s data (2003 June) corrected 
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-0.5 







0.5 



phase 



Fig. 2. — Calculated Rc light curves are plotted against the binary phase (binary phase is repeated twice from —1.0 to 1.0) together with 
the observational points of Kato et al. (2004, square) and of Kiss et al. (2004, circle). A thick solid line denotes Rc light curve for the best 
fit model. The apparent distance modulus is (m — M)^ = 12.05 for the WD luminosity of LwD = 3, 000 Lq . The adopted system parameters 
are listed in Table 1. 
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VI 494 Aql 1.0 WD + 0.3 MS (m-M)p=1 2.05 i=78.5° 

Kiss et al.'s data (2003 August) 
, Kato et al.'s data (2003 June) corrected 
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Fig. 3. — Same as in Fig. 2 but for variations of numerical models with different heights of spirals. The variations appear at orbital phase 
of 0.1 — 0.3 but the peak at phase ~ 0.6 — 0.7 is robust. See text for details. 



